Preprint WSU-NT-14-1996 



QCD evolution with longitudinal fields and heavy quarks 

A. Makhlin and E. Surdutovich 

Institute for Nuclear Theory, University of Washington, Seattle, WA 98195 

and 

Department of Physics and Astronomy, Wayne State University, Detroit, MI 48202 

(October 31, 1996) 

QCD evolution equations that naturally include longitudinal (non-propagating) fields and heavy 
quarks are derived. We start with the integral equations of quantum field kinetics and obtain the 
master equations, similar to DGLAP evolution equations after several consecutive approximations. 
We demonstrate that in their primary form, the evolution equations include a new element, feed- 
back via longitudinal fields, leading to a low-a; enhancement in the e-p DIS cross section. We show 
that the structure function Fl is very sensitive to the dynamics of the longitudinal fields and that 
the heavy quarks in evolution equations make this effect even more pronounced. 



I. INTRODUCTION 



In this paper we address the dynamics of longitudinal fields and heavy quark fields in high-energy processes, keeping 
in mind practical importance of these questions for the theory of ultra-relativistic heavy ion collisions. One of the 
main difficulties of this (not yet existing) theory is the impossibility of describing the initial state (two nuclei) , the 
transient process (which, hopefully, includes the QGP stage) and the hadrons of the final state using the standard 
language of scattering theory. In fact, even in the much simpler problem of the pp-collisions, one has to make certain 
model assumptions (like, i.g., the parton model) in order to put the problem in the form required by the scattering 
theory. Unfortunately, for nuclear collisions, the situation is much worse, since one cannot even think about consistent 
decomposition of the whole process into the hard and the soft parts; the factorization scheme is inapplicable. The 
most unfortunate circumstance is that the structure functions of the protons, the partonic substitutes for the unknown 
wave functions, cannot be used any more. 

Two possible ways to cure the problem were suggested recently. McLerran and Venugopolan have started with a 
classical model of a large nucleus in the infinite momentum frame (IMF) . This model has been gradually improved 
by Jalilian-Marian, Kovner, McLerran, Venugopolan and Weigert by accounting for the small quantum fluctuations 
in the strong "external" classical field |^||. A remarkable yield of this ongoing study is an understanding that the 
so-called low-a; enhancement may be an effect of the classical field itself Q. 

The idea of a geometrical and dynamical similarity between all deeply inelastic high energy processes, including the 
inelastic high energy ep-scattering, is the guide for our study IH. The method of Quantum Field Kinetics (QFK) 
H allows one to show that the dynamics of all inclusive processes is sequential on the real time scale and ends at the 
moment of measurement. Consequently, one can avoid factorization in the derivation of the evolution equations [0. 
The best prospect of this approach is connected with the hope of extracting from the DIS data more rich information 
(e.g., information about the dynamics of the static fields,) than is usually done with the aid of the DGLAP evolution 
equations of the parton model The causal character of the QCD evolution which is clearly visible due to the 

QFK technique does imply such an opportunity. As a result, one may treat the hadronic collision not as the scattering 
of partons, but as the interaction of the fields from the colliding hadrons or nuclei. In this paper, we continue the 
study of Ref. ||] and show that the presence of the static (not propagating) fields in the evolution equations leads to 
a steep power-like behavior of the inelastic ep-cross section at low x. 

The QFK approach to deep inelastic processes allows us to address the question of the accuracy of the evolution 
equations as well. The complexity of this problem in the calculations based on the operators product expansion (OPE) 
is well understood. Ranking the operator functions by their twists and the coefficient functions by the order of their 
perturbative expansion, the OPE employs two essentially different expansions, and their interplay is difficult to keep 
under the control (see, i.g., Ref. |lC|]). Our approach is different; we descend from the Schwinger-Dyson equations 
all the way to the DGLAP equations and study what of physical importance is lost on this way. We show that the 
dynamics of the classical (longitudinal) fields is an unavoidable partner of the propagating (transverse) fields in the 
QCD evolution. As a result, the equations for the transverse fields acquire a kind of feed-back via the inherently 



1 



inseparable longitudinal fields. By examination, this new element may alone lead to the low- a; enhancement of the 
inclusive cross-section. One does not even need to have a classical source like the valence quarks of the McLerran- 
Venugopalan model. In order to obtain the master equations, like DGLAP, one has to exclude the longitudinal fields 
from the evolution equations by the brute force. We argue that this step is equivalent to the introduction of the 
parton picture. It may be consistent only if the so-called factorization scale, /ip, is a measured parameter. 

In their full form, the evolution equations are nonlocal in the transverse direction. In conjunction with the collinear 
problems of null-plane dynamics, this leads to severe infrared singularities, which are easily regulated physically, but 
at the price of abandoning the null-plane dynamics as a technical tool. In order to obtain the master equations, any 
evidence of the non-locality has to be eliminated, once again, by a brute force. Thus, we solve (or at least address) 
the problem of the accuracy of the DGLAP equations by the extension of the system of evolution equations beyond 
their standard mathematical form and by a revision of their physical content. 

The way we derive and study the evolution equations allows us to incorporate the massive quark fields into the 
QCD evolution. This is a very real problem since a vast set of the current and forthcoming experiments requires 
a theoretical understanding of the role of heavy quarks in high energy processes. At HERA energies, the range of 
the accessible momentum transfer between the electron and the proton is very wide and the c- and 6-mesons are 
quite frequently met among the secondaries, affecting the measured cross section of deep inelastic scattering (DIS). 
A significant multiplicity of charmed quarks is expected in heavy-ion collisions at RHIC energies. Their evolution 
may carry important signatures if the QGP was a part of the entire scenario. Here, the theoretical calculations rely 
heavily on the structure functions which are taken from DIS data. The definition of the heavy flavor DIS structure 
functions (or sources) should be identical to those of gluons and light sea quarks. Only in this way can one consistently 
introduce the notion of intrinsic charm and beauty, and describe their excitations unambiguously. The contribution 
of heavy flavors to the ep-DIS cross section should naturally die out when we move in the direction of lower Q 2 along 
the evolution scale. However, the status of heavy quarks in QCD evolution equations has remained uncertain for over 
two decades. To handle this problem, one may wish to rely on a minimal extension of the renormalization group 
method and account for the heavy quarks only inside the fermion loop of the gluon self-energy, as is done in Ref. 
|pr| . A kind of transient regime in the vicinity of the partonic threshold Q 2 ~ m 2 was suggested in Refs. |l2l , |l3f . In 
these calculations, various regions (below and above the threshold) were treated in a different manner: with different 
numbers of active flavors, and different renormalization prescriptions. Special subtraction schemes were introduced at 
the threshold. However, this method does not eliminate the problem itself: special separate treatment of the intrinsic 
charm and beauty, the wee parton inside the proton, and the extrinsic charm and beauty, the heavy quarks created 
in 7*g-interaction. As it will be seen later, both these treatments of the massive quark fields in QCD evolution are 
incomplete and it is expedient to consider this problem together with the problem of the evolution of static fields. 

II. DEEP INELASTIC ELECTRON-PROTON SCATTERING 

Our approach to high-energy collisions amounts to the elaboration of a way to incorporate (to extract and to use) 
the dynamical information about the transient processes which takes place during the collision. Therefore, we start 
with the simplest example of inelastic ep-scattering where at least one of the participants is structureless and pay 
special attention to the procedure of measurement. We divide this section into the three parts. For the sake of 
completeness we begin with a brief definition of the DIS cross-section in terms of quantum field kinetics and define 
the null-plane variables which will be used for all the following calculations. Next, we discuss the temporal sequence 
of dynamical processes during inclusive measurements. Finally, we derive the explicit formula of the measurement 
which does not imply the parton decomposition of the proton and accounts for the static component of the proton's 
wave function. 



A. Observables for DIS 

Let us consider a collision process with the parameters of only one final-state particle explicitly measured. Let this 
particle be the electron with momentum k' and spin a' . The deep inelastic electron-proton scattering is an example 
of such an experiment. All the vectors of the final states which are accepted into the data ensemble are of the form, 
a^/(k')|X) where the vectors \X) form a full set. The initial state consists of an electron with momentum k and spin a 
and some other particle or composite system carrying quantum numbers P. Thus, the initial state vector is aJ.(k)|P). 
The inclusive transition amplitude reads as (X\a a i(k') S aJ.(k)|P) and the inclusive momentum distribution of the 
final state electron is the sum of the squared moduli of these amplitudes over the full set of non-controlled states \X). 
Therefore, we obtain the following formula, 
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^ = (P\a a (k)S^al(k')a a/ (k')Sai (k)|P) , (2.1) 

which is just an expectation value of the Heisenberg operator of the number of final state electrons over the initial 
state. Since the state |P) contains no electrons, one may commute electron creation and annihilation operators with 
the ^-matrix and its conjugate, , pulling the Fock operators a and to the right and to the left, respectively. Let 
one-particle wave function of the electron. Then the procedure results in the following formula, 

™1 = I y [ dxdx'dydy'^ix^Ux 1 )^^^ ( J S \ 3 ) \P)^ + J \y)^Uy') . (2.2) 

Introducing the Keldysh convention about the contour ordering [[s||l5|l, we may rewrite this expression as 



„T = -YY(-1) A+B f dxdx'dydy'^J(x)^l(x')(P l - 

where S c = S. The electron couples only to the electromagnetic field. Therefore, to the lowest order in electromag- 
netic coupling e, 

= yE/ dxdy^l(y)^i + J(x){P\ A(x) A(y)\P)^ + J (y)^l(x) , (2.4) 

where A(x) is the Heisenberg operator of the electromagnetic field. Already at this very early stage of the calculations, 
the answer has a clear physical interpretation. Since only the final state electron is measured, the probability of the 
electron scattering is entirely defined by the electromagnetic field produced by the rest of the system. The fina l 
answer can be obtained either by using the equations of QFK |qJj, or one may take short cut and iterate Eq. (2.4) 
by means of the Yang-Feldman equation, A(x) — J d 4 yA ret (x,y)j(y) , where j(y) is the Heisenberg operator of the 
electromagnetic current and A ret (x,y) is the retarded propagator of the photon. 

Summation over the electrons spins brings in the leptonic tensor L^ v (k, k'). If q = k — k' is the space- like momentum 
transfer, then the DIS cross section is given by 

^ = (^ h ^ 1 KtMW a ,{ q )^ , (2.5) 

where W^ v {q) is the standard Bjorken notation for the correlator of two electromagnetic currents, 

W» v (q) = ^^[-in%(q)] . (2.6) 

normalized in a way which provides correspondence to the parton model, and tt^(x, y) = — i{P\]^{x)]" (y)\P) is the 
electromagnetic polarization tensor, correlator of the currents, which are the sources of the field which has scattered 
the electron. We accept the standard decomposition of the correlator W^ v in momentum space, 

where v = qP, Q 2 = —q 2 > 0, xgj = Q 2 /2v and 

^ = -<r+*£; C ^ = -^ + ^ + ^- g ^. (2. 8 ) 
q z v v l 

We shall perform all computations using the infinite momentum frame fixed by the null-plane vector n^, = 
(l,0 t ,— 1), n 2 = 0. It defines the "+" -components of the Lorentz vectors, na = a + = 2a_ = a + a 3 ; a~ = 
2a + = a — a 3 . In the infinite momentum frame, the 4- vector of the proton's momentum has components 
pp- = (P+/2, t , P + /2), P~ = P° — P 3 = 0. The vector of the momentum transfer has the components, 
q» = (^/P+,q t ,-^/P + ), q + = 0, q~ = 2v / P+ . 

Instead of the invariant W2, we shall use the mass- independent structure function F2(xBj,Q 2 ) — VW2/M 2 , which 
is calculated via the equation 

(P+) 2 F 2 

c 2 = W^n^nu = K ' 2 . (2.9) 

The longitudinal structure function, F^XBj ,Q 2 ) = Wl, should be calculated in accordance with 

3P L = 2x BjCl + 2F 2 ; c x = W^g^. (2.10) 



(2-3) 
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B. Causality and temporal order in DIS 



The position of the two photon propagators, A ret (q) and A ae i v (q) ( since A a d v enters with inverted space-time 



arguments, both are retarded propagators ) in Eq. (2.5) reflects the two major aspects of causality in relativistic 
quantum mechanics. First, any statement concerning the time ordering is in manifest agreement with the light cone 
boundaries. Second, the dynamical information about the quantum- mechanical evolution is read out only after the 
evolution is interrupted by the measurement ( the wave function has collapsed). Therefore, the standard inclusive 
e-p DIS delivers information about quantum fluctuations associated with the proton. The space-like photon which 
scatters the electron "belongs" to the proton and is, in fact, its part. The method of quantum field kinetics H allows 
one to continue this line of reasoning by employing the basic definition of the electro mag netic polarization tensor [^) . 



Then, neglecting any corrections to the electromagnetic vertex, we may rewrite Eq. (2.6) in the following way, 

W^(q)=^e 2 f ^fP^ J J0- A TrYGU P + q)YGl 1 {p) . (2.11) 

In what follows we shall calculate only the non-singlet functions of a given flavor and omit the flavor label / when 
it causes no confusion. The quark field correlators in this equation form a 2 x 2 matrix and obey a matrix integral 
equation of the Schwinger-Dyson type, 

Gab = Gab + GarX'RsGsb , (2-12) 

which allows for a symbolic solution ||. This solution expresses the entire matrix of field correlators in terms of the 
full retarded and advanced propagators and the sources Y*ab (the "current" correlators), 

Gab = G re tG ^G abG ^G a dv + (— ^) A+B G ret T,ABG ad v > (2-13) 

where GT^. and GT^. are the differential Dirac operators acting on the right and on the left, respectively. The retarded 
and advanced Green's functions obey more familiar equations, 

G« t = Gret +G„ t E« t G„ t , (2.14) 
which allow for the symbolic solutions also, 

G;i =G-} t -E ret . (2.15) 



A simple examination of Eqs. ( 2.13 ) shows that all four dressed correlators Gab can be found as formal solution of 
the retarded Cauchy problem, with the bare correlators Gab as the initial data and the self-energies as the sources. 
Indeed, integrating the first term of each these equations twice by parts, we find for all four elements of Gab ■ 

Gi (x,y)= f d^d^G ret (x,Ol fi G 1 o(^vh v G adv (r,,y) - f d^d 4 r,G ret (x,0^,v)G a dv(v,y), (2.16) 
01 J 01 J 01 



Goo(x,y) = (d^d^Gr^x.O^Goo^YGad^y)- fd^d A r,G ret (x, £)[±<V + Eu(& v)]G a dv(v, v)- (2-17) 

Though in the original derivation of these equations we employed the "bare" fields of the interaction picture, the 
last two formulae allow one to make a step further and to replace the correlators Gab of the bare fields with the 
actual values of correlators Gab on the hypersurface of the initi al data which, in their turn, were created in the course 
of the preceding evolution. After this modification, Eqs. ( 2.15|) — ( 2.17 ) become an excellent tool for the study of all 



inclusive inelastic processes. Similar equations can be derived for the correlators of the gauge boson fields and their 
self-energies. 

Viewed in the context of the QFK, the physical picture of the QCD-evolution looks like the virtual assembling of 
the source of a proper electromagnetic fluctuation, one which is capable of providing a given momentum transfer to 
the electron. In the event such a transfer is detected, the hadronic system has to recoil and, provided the transferred 
momentum is high enough, to emit a backward quark jet. This jet is a successor of the free quark in the perturbative 
vacuum. It is a general principle of the quantum mechanics of radiation that the emission is the excitation of the 
field mode. To be excited, the mode itself should exist. The normal modes of QCD are hadrons and they have a true 
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physical nonperturbative vacuum as the ground state. The process of the emission of a free quark is impossible without 
creation a new ground state. In the real world of asymptotically stable states, the words "perturbative vacuum" mean 
a special highly excited state of the hadronic system. The QCD evolution includes the description of its creation by 
definition. Thus, we are unavoidably led to a physical picture based on the long standing observation by Shuryak |lq | 
that the energy scale of confinement is weak, while the energy density of the QCD vacuum is very high. From this 
point of view, multiparticle production in hadronic collisions looks like a signature of a "cavitation" of the physical 
QCD vacuum and the creation of a domain with a perturbative dynamics of quarks and gluons. To switch it on, one 
needs extreme Lorentz contraction which resolve fluctuations of much shorter scale than those providing confinement 
in hadrons. The entire process is strongly localized in space and time ]Ti] ). To describe this process quantitatively, we 
have to include two types of correlations: 

Gio = G 10 — G re tYiioG a dv j Dm = Df" a — D re tHioD a d v , (2.18) 

01 01 01 01 01 01 

where the superscript "#" is assigned to all the states in the continuum of the free on-mass-shell fields, 

Gfo\p) = -2TTi5 l:j (0 + m)9(±p o )S(p 2 -to 2 ), (2.19) 

01 

D* ab '^(p) = -2mS ab d^(p)0(±p o )5(p 2 ) , d^(p) = -<T + V - i!LE_ . (2.20) 



(Projector d^ v {p) corresponds to the null-plane gauge A + = 0.) These states are initially empty and the vacuum 
correlators Gf 001 (p) and Df a Q1 (p) represent only on- mass-shell quarks and gluons in the final states. The second 
terms in Eqs. (|2.18|) include the off-diagonal self-energies Si and flio, the sources, which are the remainders of 

' 01 01 

the classical configuration of fields in the hadron and in the physical vacuum at some moment of the dynamical 
evolution before the final interaction (the measurement). Before this interaction happens, the decomposition of the 
field configuration into the radiation field and the classical remainders is virtual; the phases of the partial waves are 
balanced to form a hadronic correlator propagating through the condensates of the physical vacuum. 



C. Kinematics of the inclusive DIS measurement 



We do not employ the method of operator product expansion where XBj immediately appears as the natural variable 
corresponding to the light-cone momentum. Moreover, we are going beyond the parton model and intend to include 
new dynamical elements, e.g., the static fields, in the evolution equations. Therefore, we have to pay special attention 
to the formula of measurement. In particular, we have to establish a connection between the kinematic variable XBj 
and the Feynman variable xp without any reference to the parton model. The ep-DIS measurements rely on a set of 
rare events where the space-like photon 7* is created by a strong change of the state of one quark only. This quark 
almost instantaneously decouples from the bulk of the proton and initiates the jet. Thus, we have to require that one 



on-mass-shell quark appears in the final state in the perturbative vacuum, and, therefore, Eq. (2.11) can be rewritten 

as 

W^(q) ^Y, ^ 1 ^ J (2^ Tr ^ ( ^ + d + m fh" G ret(p)Zoi(p)G adv (p)] [2mS[(p + q) 2 -to 2 ]] . (2.21) 
By virtue of Eq. (2.9), and introducing x = xf = P + / P + , we obtain the expression, 

F 2 (x Bj ,Q 2 ) = e /^^F I xdx J d 2 Ptdp-6[p+p- + 2ux - (p t - q t ) 2 )] ip+T^{p) , (2.22) 

where J 7 ^' is the coefficient in the following decomposition of the quark source, 

G rei £ 10 G acfo = |(?<S)(p)+fljj + flP(p)+mCP(p)+wP(p)(|J^-rijO . (2.23) 
01 

Following Ref. we can express the invariant J-{p) via the "unintegrated structure function," 
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where v = Pq. [As it will be seen later, in a certain approximation, q/(x, Q 2 ) becomes a structure function of the 
deep inelastic scattering of a- quark of flav or /.] Introduction of the (5-function 8(p~) at this step of calculations is 



motivated physically. According to Eq. (2.5), the field fragment with the momentum p M belongs to the yet undestroyed 
proton. Let f{x + ,x~) = f{x° + x 3 ,x° — x 3 ) be some proton-related field quantity, e.g., the proper electromagnetic 
field of the proton. If the proton, before its interaction with the electron, is viewed as a kind of a wave packet that 
propagates without dispersion at the speed of light, then the function / should depend only on a single argument x~ . 
In the momentum representation, we have 



f(x + ,x )= / dp+dp f(p + ,p )e = l 

Thus, in order to eliminate the dependence on x + and to keep the proton as a stable wave packet which propagates 
along the light cone, one should require that f(p + ,p~) = f(p + )8(p~)- Admitting the opposite, we would run into a 
conflict with the Lorentz contraction of the initial state proton, viz., the radiative corrections to its "valence structure" 
would increase its size. The dynamics of the quark and gluon fields behind this condition is totally nonperturbative. 
In the same way, following Eq. ( p. 10 ), we obtain the expression for the longitudinal structure function, 



ie 2 \f f 

3F L (x Bj ,Q 2 ) = / dx / d 2 p t dp-8[p + p- + 2vx - {fit - q t ) 2 )]Fi(p) , (2.25) 

n ( 27r J Jo J 

where 

T\ = [2v(x% - x 2 Bj ) - x Bj {p 2 + m 2 )]T - 2{p+) 2 x Bj B + Am 2 x B] C. 

The next calculations are as follows. We integrate over the angle between the vectors fit and qt and re-scale p\ in the 
remaining integrals by 2v, y — p 2 /2v. Denoting \i — m 2 j2v , we obtain, 

F **»^ = T f XdX f^^'dy ^f^- ]2 nl/2 • (2.26) 



Here, the denominator is the quadrupled area of the triangle with the sides pt, qt, and V2 



vx 



by 2v. The limits of integration in the Eq. (2.26) come from the triangle inequality. Now, we should take the limit 
v — * oo keeping the x B j finite. In the new variables, the parameter v ent ers o nly the argument of the function 
dqj{x,2vy) / dy. This function vanishes in the limit of v — > oo. The integral ( [2.26| ) vanishes also unless the variable 
y is allowed to take arbitrary small values. Therefore the lower limit of the integration over y should be set equal to 
zero, i.e., 

x F = x Bj +A* = x B j(l + m 2 /Q 2 ) . (2.27) 

For massless quarks, we obtain the usual condition, xf = x B j. [ If we were doing the same analysis in the natural 
momentum variables, then the lower limit of the integration over pt becomes infinite when v — > oo at fixed x B j. In 
fact, the finite interv al of integration over p t would keep its length 2Q and slide as a whole to infinit ely h igh values of 



p t . However, if Eq. ( 2.27 ) is satisfied, then both limits become independent of v. ] The condition ( 2.27 ) defines the 



upper limit of the p t -integration integration also, and the final result reads as 

F 2 (x Bj , Q 2 ) = e) (x Bj + f 2 dp 2 dg/(**i+yV2*,P?) (2 . 2g) 

Here, the shift of x B j by m 2 /2v is a standard quark mass correction Change of the upper limit from the 

familiar Q 2 to AQ 2 is unusual but not unexpected. Indeed, of the three sides of the triangle, p t , qt and \/2vx — m 2 , 



the third one, by virtue of Eq.(2.27), is equal to qt. Therefore, the length of the vector p t , the transverse momentum 
of the incoming off-mass-shell field, varies from (at p t — qt) to 2q t (at fit — —qt). Furthermore, with logarithmic 
accuracy there is no theoretical difference between log Q 2 and log AQ 2 when Q 2 — > oo and no conflict with the leading- 
log approximation (LLA), the AP-equations, is anticipated. [For the usual on-mass-shell parton, we put p t — and 
p~ = and obtain Eq. ( p.27| ) directly from the delta-function in the integrand of Eq. ( 2.22[ ) . 



As it was discussed in Ref. |6[, the derivatives of the structure functions like dq{x, Q )/dQ 2 (the sources), are 
similar to the densities of states, p(E), of the condensed matter or scattering theories, and the structure functions 
themselves are similar to the number of states n(E) below the energy E. In an experiment we measure n(E) 
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(which is proportional to the kinematically allowed volume in the phase space), rather than p(E). Therefore, it is 
quite natural that the upper limit of the pt integration coincides with the kinematic boundary 2Q. The less trivial 
fact is that the proof of the resonant condition of the m easu rement, xbj — xp, imposes a limit on the rate of the 
QCD evolution at h igh Q 2 ; the integral over p\ in Eq. ( 2.22 ) should be convergent. One more important aspect of 
equations like (2.28 ) is that the differential dqj(x,p 2 ) works like a measure in the functional space of the of the 
quantum- mechanical fluctuations of the quark field. Below, a similar kind of measure, dG(x,p 2 ), will be introduced 
for the gluon field. However, the measure will never appear for those components of the fields which do not have the 
property of propagation and are not independent variables in the Hamiltonian formulation of the system's dynamics. 
These are the static components of the fermion field and longi tudin al gluo n fiel d. 

The full expressions for the invariants T and T\ in Eqs. ( 2.22 ) and ( 2.25 ) are very long and their exact form 
exce eds th e accuracy of approximations used in this paper. The complexity comes from the product G rei EoiG a( fo in 
Eq. (2.21). Indeed, we have two vectors, p 11 and n 11 , at our disposal. The fermion field is massive and none of the 
fields are polarized. Thus, all the self-energies ~Eab have three terms in their spinor decomposition, 



S(p) 



m)a 2 (p) + jhp + a 3 (p) + ma m (p) 



(2.29) 



The scalar invariant functions of this decomposition can be easily found from the three Dirac traces, 
TrE = 4m(a m (p) - a 2 {p)), Tr[^E] = 4p + <J 2 {p), Tr[pT] = 4[p 2 a 2 (p) + (p+) 2 <7 3 (p)] . 



(2.30) 



By virtue of Eqs. (2.15) and (2.29), the retarded and advanced propagators are readily found as 



Gre 

ad 

where 



[{$-m){l-a\ A ')- iP o \ 



a), 



ftp a 



(2.31) 



Q (l) 



(p 2 -m 2 -aP)(l-aP) 



m 2 \aP] 2 



[sP] 



(2.32) 



and the following short-hand notation is introduced: ctq = (p — m )fa + 2(p ) 03 + 2m a m , 



; (2) 



(p) = p 2 ~ m2 



aP (p) , sP (p) = 1 - Jp (p) , sP (p) = maP (p) 



For the massless quark field, we can proceed without approximations and obtain a reasonably simple expression, 



G ret (p) E ^01^ (p) Gadv (p) 



[/p- MP 2 /2p+)}ap°\p) , jh ap ] 



(P) 



S${p)S£{p) 



2p+ S?(p)Sf(p) ' 



( m = ) 



(2.33) 



which explicitly exhibits separation between the dynamical and th e con strained (x + -instantaneous or static) parts of 
the fermion field. The constrained part (the second term in Eq. ( [2.33] ) ) is easily recognized since (after the string 
G re tTiGadv is assembled) it has no pole corresponding to the propagation. A reasonable approximation for the massive 
quark field will be to neglect the radiative corrections in the retarded and advanced propagators. Then 



T u \p) 



If) 



IP 2 



,212 



B (f) 



,(/) 



C U) 



If) 



If) 



(2.34) 



and the expressions for the functions F 2 and Fl become significantly simpler. The factor [p 2 — m 2 ]~ 2 in the above 
formulae is never singular; it originates from the product of the retarded and advanced Green functions and should be 
read with the principal value prescription. With this simplified form of the invariants, the integrand of the structure 
function F 2 , which is defined by Eq. J2.28D , takes the following form, 



dg f (x,p 2 ) _ V lab P+ , a 

"yP > J„2 — n/r,_\S l P 



i f) (p) 



dp 2 t 



2(2vr) 3 [p 2 ~ m 2 } 2 ' 
The integrand of the longitudinal structure function Fl becomes relatively simple also, 

,2 Jf), 



2m) cjtf'ip) 



2m 2 f 



[p 2 — m A 



212 



p 2 -m 2 ° m 



U) (p)~4 f \p) , 



(2.35) 



(2.36) 
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and the final expression can be conveniently written down in the following form, 



„ , 9x 2to 2 Q 2 (4Q 2 + to 2 ) 

3Q 2 F L (x Bj ,Q 2 ) = f Jn2 \ fJ F 2 (x Bj ,Q 2 ) 



[Q 



2e 2 f x Bj Q 2 r*Q 



+ m 



f J o 



d P 2 t 



dM f (x Bj (l + m 2 /Q 2 ),p 2 ) 



dp 2 t 



+ Q f (x Bj (l + m 2 f /Q 2 ),p 2 ) 



where we have introduced the following notation, 



2(2tt) 3 " p 2 t -m 2 



2(2tt) 3 



(2.37) 



(2.38) 



The first term in Eq. (2.37) is in direct proportion to F%. It is entirely due to the finite quark mass to/. It shows that 
at Q 2 r~j m 2 the non-singlet structure functions F 2 and Fl of the heavy quark may be of the same order. 



III. EVOLUTION OF THE SOURCES 



Our next step is to find evolution equations for the sources of the electromagnetic fluctuation which scatters the 
electron. These sources are given by the off-diagonal self-energies Eoi (for the quarks) and Eio (for the anti-quarks). 
Their general expressions were derived in Refs. [pl,pfl and, as in Rcf. §, we shall use approximation of the bare vertices. 
Furthermore, we shall use simplified equations ( |2 . 1 8J ) which define evolution of the quark sources and make similar 
simplifications for the gluons (hereafter, the gauge A + = for the gluon field is assumed). 

The general expression for the gluon sources, Ilvon, were derived in Refs. [^|^] also and an approximation of the 

bare vertices will be used here as well, 



d A k 



E(or)(p) = ig r C F J ^y I {7 A1 G re t(fc)E(-oi)(fc)G ad ,,(fc)7 ! ,I)^oj I (fc - p) 



(3.1) 



DoiGO = ^9 2 {^ J 0^Tv{^G ret (k)J: Q1 (k)G adv (k)YG*(k- p) + ^G*(k+p) 7 u G ret {k)E 01 (k)G adv (k)} 

+ \ J ^V a ^{p,k- P ,k)[D ret (k)n 01 (k)D ad 4k)] a J^ (3.2) 

The renormalization of these equations as well as the study of their infra-red behavior is a special problem. The idea 
of renormalization in the QFK-based calculations has been discussed in Ref. [f| for the case of pure glue-dynamics. In 
this simplest case the coupling constant should be replaced by the "running coupling," g 2 — > 4ira s (p 2 ). However, we 
still have not obtained satisfactory balance of the collinear singularities except for the limiting case when the evolution 
is described by the DGLAP equations. The status of this problem will be discussed separately Q. 



A. Splitting of the evolution equations 



The polarization tensor IF 1 " only appears between the retarded and advanced propagators, i.e., in the combination 
[D ret (p)I\-(p)D a dv{p)\ tJ ' v ■ Both propagators contain the same projector, d^(p) which (by the gauge condition) is 
orthogonal to the 4- vector n M . So, of the general tensor form, only two terms survive, 

rF» = gT p 2 Wl (p) + p^w 2 (p). (3.3) 

The others, like + n^p v or n^n v , will cancel out. Introducing one more projector, 
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which is orthogonal to both vectors n v and p^ , we find that the invariants Wi and w 2 can be found from two 
convolutions, 



- d ia ,{p)W v {p) = 2p 2 Wl (p), and nil n v W v (j)) = [p+f w 2 {p) , 



(3.5) 



independently of the other invariants accompanying the missing tensor structures. The new projector, which includes 
only two transversal gluon modes, naturally appears in the tensor with a gluon source. Indeed, solution of the 
Schwinger-Dyson equation for the retarded gluon propagator, 



D„ t =D ret + D rct U r<st T> rct . 



(3.6) 



can be cast in the form, 



D S:.) <rt 

(R 



p 2 -p 2 w { p(p) ( p+ ^ l-wP(p) 



(3.7) 



(2) ( 



.(2) 



With the shorthand notation, Wf (p) = p —p w\ (p), and W 2 (p) = 1 — w 2 (p), we easily obtain, 



D ret (k)ILfoi\ {k)D a dv(k) 



-df" \p)p 2 w ( p (p) 4 ? ° ) 



(p)r 



w^CpWCp) (p*) 2 w 2 fl (p)vi#(p)' 



(3.8) 



Once again, similar to the case of the fermion field, we encounter the longitudinal (x + -instantaneous) part w 2 of the 
gluon source which is not accompanied by the propagation poles in the string D re JID a dv 

The role of the tree propagators and which include the radiative corrections is clearly understood. These 

corrections keep memory about the balance of the phases between various fields in the yet undestroyed proton. In 
fact, they nonperturbatively maintain the <5(p _ )-prescription, the requirement that the proton does not fall apart in 
the course of the QCD evolution. Keeping this important qualitative observation in mind we shall limit ourselves, in 
what follows, to the bare retarded and advanced (tree) propagators of the quark and gluon fields. 

Our next goa l is t o find the evolution equations for the three invariants, a^\x,p 2 ), (x,Pt) an d ffm (%,Pt)- 
Using the Eqs . ( |2.30) an d (3.5) we easily extract evolution equations for various invariants from the tensor evolution 
equations (|3.l|) and (|3.2[). In sequence, we obtain three equations for three invariants of the fermion source, 



4°w = W) 







p — mj 



z 2 + l 



- 2m 



[k 2 - 



(2tt) 



d 4 k\ C F S+[{k-p) 2 

<j { i\k) ■ 



,212 



+ za 2 J '(k)-2mj(l- z) 



k 2 



+S+[(k - p) 2 - m 2 } 



(/)^_gg!M) 



[(• 



TYl j 



(l-z) 2 + z 2 



,(io) 



+ m 2 j\ k 



2 ^1 



(-*) 



[k 2 



a^(p) 



(27T) 



d 4 k{ C F 8 + [{k-p) 2 } - 



a { /\k) 



z [k 2 



m 



212 



+ 8+[(k-p) 2 -mj 



(27T) 



d 4 k\ C F 8+[(k-p) 2 



z(l- z)w { 2 W) {-k) 

l-z k 2 wt 0) {-k) 
2z 

<T { i ] {k) 



(3.9) 
(3.10) 



[k 2 



m 



212 



/J ~ m f 
1 fc 2 ^ 10) (-fc) 



+S+[(k - pf - m 2 } 



2z 



[k 2 



(3.11) 



and two equations for the invariants of the gluon source, 



v 2 w m (p)= 2g2( ^ } 

P W l W (27r)3 



-N c 5 + [(k-p) 2 } 



mj) 



(z-l) 2 + l 



2m 



J d^Tf^S+Kk-p) 2 -™ 2 
{f) (k) 



z 

1 - z 

z 



fi [k 2 



n 2l2 



+ (l-z)a 2 J> (k)-2m}z-^ 



1 - z 



+ z(l-z))k 



(01) 



^- (z--) 2 w (01) (k) 
[k 2 ] 2 [Z 2 )W2 (k) _ 



(3.12) 
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w. 



(or 



(P) 



(2tt) 3 



f 



,(/) 



(k) 



212 



(3.13) 



Here, the sum ov er (/ ) runs over all quark flavors and anti-flavors. The longitudinal fields appear in the evolution 
equations (2.13)— (2.18) in an alternating regime. This remarkable feature of the evolution equations has a very clear 
physical explanation. Indeed, if the static source has emitted a propagating wave, at least one additional emission is 
necessary to create a new stati c fiel d config uration. 

Deriving the field equations ( 2. 13| )— ( 2. 18 ) we were relying upon the most straightforward implementation of the idea 
of quantum mechanical evolution; viz., the Heisenberg picture of the evolution of obse rv ables of DIS with the quark 
and gluon fields of a stable hadron in physical vacuum as the initial data. Eqs. ( |3.9| )-( 3~13| ) are the mathematical 
expression of this physical picture. The initial data are not given explicitly, and if they were, the range of the 
incorporated physical information should be the same as in the local operators of the OPE. However, we do not view 
the perturbative part of the evolution as the scattering problem and this results in the new feature of the above 
equations; the invariants a 2 and w 2 which correspond to the sources of longitudinal fields participate in the evolution 
equations on the same footing as the sources cto and w\ of the transverse fields. Therefore, we have a reason to 
suspect that the entire dynamics of the QCD evolution carries a classical component. This is in compliance with the 
qualitative understanding of the QCD evolution as the virtual decomposition of the hadrons, which are the genuine 
fundamental modes of QCD, in terms of the alien set of modes defined as excitations above an artificial perturbative 
vacuum. The parton picture should emerge after certain approximations, different from the formal twist classification 
of the local operators of the OPE expansion. 



B. Feed-back via longitudinal fields 



In what follows, we shall make several approximations and try to estimate their accuracy. In the OPE-based 
calculations, it is a formidable task. Indeed, the hierarchy of scales in OPE is two-fold. On the one hand, the operator 
functions are ordered by their twist which, would we wish to treat it statistically, corresponds to the order of the 
irreducible correlation function in condensed matter physics. On the other hand, the coefficient functions are ranked 
by the order of their expansion in powers of the coupling constant. Because of the capricious interplay of these two 
essentially different expansions it is difficult to obtain a reliable estimate of the accuracy even for the widely used 
DGLAP equations. 

We shall pose the problem in a differ ent way and try to understand what elements of the physical picture are 
lost when the field equations (3.9)— (3.13) are approximated by the master equations of the parton evolution. Let us 
consider the transverse components of the sources, <jq and w\, as granted and use the connection to the observables, 
equations ( J2.36j )— (p73^) , 



c ip 



[p 2 



6(p- 



dp 



cip+4 f) (p)=5(p-)Q f (x,p 2 t ) , 



for the fermion field, and introduce the similar links for the gluon field, 



c ip 



+ p 2 wi(p) _ _^dG(x,pf) . + _ f ^ _ z,„-\ n , m ji\ 



[P 



212 



S(p-)- 



dpt 



, c ip + w 2 {p) = 6(p )G(x,p t ) , 



(3.14) 



(3.15) 



where c is an insignificant common normalization constant. Integrating Eqs. ( |3 .10 ) and ( |3.13 ) over p and using the 
delta-functions to integrate the variable k~ out, we obtain, 



Qf(x,p 2 t )=c J dp- ip+4%) = 



c, 



f P+ dk+ 1 f p+ dib+ , , 



(3.16) 



where 



dp ip + w 2 (p) 



9 2 (Pt 
(2tt) 3 



q f (k+)-N c 



P+ 



dk+ (2p+ - k+f 



k+ -p-t 



Ap+W 



G{k^ 



(3.17) 
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dp, 



dpi 



and G(p + ) = J dp, 



dG(x,p 2 ) 
dpi 



(3.18) 



are the £- fractions of the quarks and the glue converted into the radiation, integrated over all transverse momenta. 
The requirements for the c onver gence of these integrals are exactly the same as in the proof of the resonant condition 
for the measurement, Eq. (2.27). The physical meaning of these equations can be uncovered by examining , e.g., the 
firsty term in Eq. (3.17), the rate of the depleting of the initial reservoir of static glue is propo rtiona l to t he to tal 
number of quarks radiated above a given value Xf- The collinear divergence at k + = p + in Eqs. ( 3.16 ) and ( 3.17 ) is 
not due to the spurious pole of the gluon propagator. It appears as a consequence of the singularity of the integration 
dp~ of the isolated mass shell delta function, 8[{k + — p + )(k~ —p~) — {kt—pt) 2 ], at the point fc + = p + . This singularity 
appears only for the longitudinal modes. It is shielded when the delta-function is accompanied by the propagators of 
the transverse field with the off-mass-shell space-like momentum, i.e., when the field is slowed down by its virtuality. 
In coordinate space, this singularity corresponds to the integration of the x + -independent function in the infinite 
limits, and it is unavoidable when the wave packet is moving without dispersion in the light-like direction, thus 
depending only on x~ . To shield this singularity, one has either to decelerate the proton (to take P + ~ y/s finite), or 
to establish a lower limit of resolution for the transverse momentum. In the second case, the natural scale of hadronic 
confinement, Aqcd, has to be taken as the cut-off. Both cut-offs play the same role; they allow for the geometrical 
separation of the fields between the near and far zones. Physically, in QCD one cannot separate two waves if their 
momenta differ by less than the width of the first peak of diffraction on an object with the size ~ Aq CD . ^ n an y 
case, integration near the lower limit will invoke a large logarithm, log (s/ Aq CD ), which will over-compensate for 
the smallness of the coupling constant. However, we must realize that this singularity is a direct consequence of the 
unphysical geometry of the fields in the infinite momentum frame; the cut-off for P + , though quite understandable 
physically, is alien to the mathematical formulation which employs a fraction xf of the light-cone momentum as the 
main variable. 

More accurate estimate of the cut-off in the splitting kernels can be obtained if the process is viewed in the scope of 
the wedge form of dynamics |L7|,[L9 1 . This dynamics employs the gauge condition A T — and it has the gauge A + = 
as its limit in vicinity of the null-p ane x~ — 0, when r — > and i] — > — oo. Here, r and rj parameterize the space-time 
coordinates in fz-plane, £ = tcos1i77, z = Tsmhr/. 

For our immediate purpose, it is enough to consider only one component of the on-mass-shell gluon correlator, 



dk t d9 4sinh?7i sinh^^e 



-ik(xi -X2) 



2(2tt) 3 {k+e-ni +k~e^)(k+e-^ +k~e^) J 2(2tt) 3 k+ + {e 2 ^ + e 2 ^){k 2 /k+) 



dk f de 



h.— g— ik(xi— x 2 ) 



(3.19) 



where k + — kte e , fc~ = kte~ e and k + k~ = k 2 . Now, it is easy to see that the pole l/k + is shielded. The cut-off is de- 
fined by the location of the interaction region; for the lowest k + one should take 771.2 ~ —Y/2 « — \n{2E max / mhadr) ~ 
\n(y/s/ Aqcd)- The invariant energy of the collision, y/s, or its equivalent, the full width Y of the hadronic rapidity 
plateau, is the measure of the resolution in longitudinal direction of the interaction which initiate the deeply inelastic 
process. Any extension of the theory beyond the parton model will include them explici tly. 



The next step is to examine the feed-back which app ears when <ji and wi from Eqs. ( 3.16 ) and ( 3.17 ) are inserted 
into the right hand side of the evolution equations (3.9) and (3.12) for the propagating (transverse) components of the 
quark and gluon fields. Since these terms in the integrand depend on k t only via the coupling constant, all integrations 
except the one over k + can be carried out, viz., 



q ^q = J dp dk dk- 
ng 2 (p 2 )k+ 1 



-^-^S + [(k-p) 2 ]6(k-)g?(kt) 



[P 2 



g 2 ((k+p) 2 )dk 



[zk 2 + (1 - z){p 2 + m 2 )Y 



pi 



m 2 + (1 — z)p: 



*g 2 (p 2 )k + 
p + p 2 t 



Kg 2 (pf)k+ 



*9 2 {p 2 t)k + 
p + (p 2 + m 2 ) 
1-z 



pi 



(1 



z)p 2 t 



(3.20) 



for the "feed-back" for the propagating fields via the longitudinal modes from quarks to quarks, quarks to glue, glue 
to glue and glue to quarks , respectively. At m = 0, all these terms yield a simple and remarkable function, 



9r(Pt) 

P + P 2 t 



k+ f(k+) dk + , 



which contains a universal Weizsacker- Williams denominator corresponding to the static field of the ultra-relativistic 
charge. This induced static source is given, in its turn, by the integral over the light-cone momenta of the transverse 
fields above the currently probed momentum p + . Almost mysteriously, this dependence appears for the fermion field 
also; it corresponds to the static pattern of the fermion wave function! 
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The intensities of the induced static sources are almost independent of the transverse momenta. It may be tempting 
to treat these sources as the equivalents of the valence quarks and gluons. However, there is a major difference. 
Evolution of the valence partons has to begin through splitting kernels, and this is the key point in derivation of the 
master AP equation; while in the field approach advocated here, the longitudinal modes are coupled to the transverse 
ones in a special way which regenerate the distribution of the field inherent to the ultra-relativistic classical source! 
Moreover, in this approach, the evolution ladder is built from the above and only from the above, solely from the 
requirement of the given momentum transfer in the measurement and in agreement with the causal picture of the 
inclusive measurement. This points us to the most important observation that the classical field pattern is inherent 
to the light front evolution, regardless what were the initial data. 

The feed-back fragment of the field evolution equations can be cast in the following form, 



dp 2 t 



/■'»■ 



(2tt) 3 



dk+ 



C, 



Pt 



Qf(k + ) 



z(l-z) 2 
m 2 f + (l- z)p\ 



(3.21) 



f dG(x,p 2 t ) 
I dp? 



f.b. 



Kg 2 {Pt 

(2tt) 3 



dk+ 
~k+~ 



-N, 



0?-l/2) 2 
P? 



g{k+) - T f 



{l-zf 



znr 



+ (1 - z)p 



S/(fc+) 



(3.22) 



where the subscript f.b. stands for the additional contribution to the rate of evolution of the transverse fields from 
the feed-back via the longitudinal ones. Q and Q stand for the left side of Eqs. (3.16) and ( |3.17 ) which connect 
longitudinal sources to the intensity of the radiated transverse fields. This contribution can be estimated in the 
following way. [From here to the end of this section we shall neglect the quark masses.] 

The data, parameterized via DGLAP equations, clearly indicate that the quantity xq(x, Q 2 ) for the sea quarks 
experiences a certain growth at small x. Let us start with the most conservativ e pro be fun ction , qf(x) = k/x, and 
calculate the feed-back along the string qf(x) — > G(x) — > G(x). Using Eqs. ( 3.16 ) and ( |3.2l| ) we obtain the first 
estimate, 



p 2 dG{x lP 2 ) 



9 2 (p 2 t) 



dp 2 



qf^Q^G 



f.b. 



^T f N c g 2 f 


l dy 


1 X x 2 ' 


f 1 du K 


(2tt)6 J x 


y 


4 y y2 _ 


Jy y 


U 


i7rT f N c g 2 K 
(2tt) 6 


' l l 

12x 


Inl-ii- 

x 9 x 


f 0(1) 





(3.23) 



which clearly indicates that our probe function is bad. The feed-back response, ~ (l/ 2 -) hi(l/x), is much bigger than 
the probe signal, ~ \/x ! This is exactly what happens in the calculations of the first order quantum correction 
to the Weizsacker- Williams field . Let us take another probe function, one which exhibits more realistic behavior, 
e.g., qf(x) = kx -1 ~ x . Repeating the same calculations, we obtain, 



o 
Pt 



9 2 {P 



dG(x,p 2 ) 
dp 2 t 



irT f N c g 2 



X 2 + \ 



f.b. 



(2tt) 6 A(A + l)(A + 2)(A + 3) x^ x 



+ 0(x~ 1 ) 



(3.24) 



Thus, our second probe is very good. The low- a: dependence 



-l-A 



appears to be an eigen-function of the feed-back 



via the longitudinal fields! To confirm this statement, let us examine the feed-back along the string G(x) — > G{x) 



G(x), taking G(x) = 71 



-l-A 



o 
Pt 



9 2 (p 2 t) 



dG(x,p 2 ) 
dp 2 t 



as the probe function: 



f X dy 


'1 


X 




2 ■ 

X 




L y 


4 " 


y 


2 

y z _ 


I 



du 

.,1+A 



f.b. ( 27 

Retaining the terms which are most singular at low x, we obtain the estimate 

G^g^G 



1 

Ay 



1 

u 



Pt 



dG(x,p 2 ) 



g 2 {p 2 ) dp 



«N 2 g 2 
(2tt) 6 



f.b. 



1 -3A 
A(A + 1) 



m To 

A 2 



A^ + A + 2 



1 



4(u - y) 



16(A+l)(A + 2)(A + 3) 



(3.25) 



(3.26) 



which supports the self-similarity of the power-like enhancement at low x with respect to the feed-back. This depen- 
dence holds both for quarks and gluons. If ln(s/A 2 ) is considered to be paramctrically large, then the leading terms 
are as follows, 
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9 
Pt 



dq(x,Pt] 



9 2 {p t ) dp 



f.b. 



*9 2 (p 2 ) 1 
(2tt) 6 x!+ a 



111 



A 2 



r 2 



Nr 



4(A + 2)(A + 3) 



(3.27) 



Pt dG(x,p 2 ) 



9 2 (Pt) dp 2 



*g 2 (Pt) i 



f.b. 



(27r) e 



In — 

^l+A A 2 



C F T f n f 



TV 2 (A 2 + A + 2) 



(A + l)(A + 2) 4(A+l)(A + 2)(A + 3) 



(3.28) 



The feed-back is positive along the string G — > — > G and negative along the strings involving quarks, like 
q — > C? — > G, g — > Q — > G, etc.. The gluons tend to boost the rate of their own evolution, while the quarks slow 
the gluon evolution down. Evolution of quarks is always boosted by the feed-back. Though Eqs. ( |3.26 )-(3.28) contain 
an infra-red regulator, lns/A 2 , which is inconsistent both with the philosophy and the technical design of the 
infinite momentum frame, it looks as though the potential source of the low- a; enhancement in deep inelastic processes 
is found correctly. The QCD evolution of observables in high-energy inclusive processes includes the evolution of 
classical (longitudinal) quark and gluon fields which can be responsible for the enhancement. The power-like behavior 
at low x has been predicted long ago by Balitsky, Fadin, Kuraev, and Lipatov from the solution of the BFKL equation 
fl4| and the exponent A was explicitly found in the case of pure glue-dynamics. Considerable work is needed to find 
this exponent in our approach. Furthermore, it is not yet clear if these two approaches rely on exactly the same 
physical input. 

Our conclusion about the power-like enhancement at low x is very close to that of McLerran et al. || but it is 
motivated in a different way. First, we neither employ nor even need the valence quarks as the classical source of the 
gluon field. In fact, we keep in mind that the strongest gluon fields are due to the vacuum condensates which exist 
even in the absence of hadrons. In the state of confinement, neither hadrons nor nuclei drag the glue; they propagate 
through it. Second, the coupling between the static and the propagating fields only weakly depends on transverse 
momentum, only via a s (p 2 ). Stronger dependence might come from the propagators, but the longitudinal modes do 
not propagate. The dynamics of the measurement always leads to induced static fields which have a steep behavior 
at low x but are integrated over all pt . For this reason, accurate measurement of the longitudinal structure function 
Fl, which is directly connected to the static component of the quark field, is of extreme importance. 

The classical pattern of the light-front evolution poses a severe problem. It is well known that energy-momentum 
conservation can be obtained either in quantum or in classical theory of radiation. Quantum theory in the presence 
of classical external field does not allow for the consistent formulation of the momentum conservation since even the 
transverse classical field can be presented at most as a superposition of states with various numbers of quanta in every 
mode (coherent states). Longitudinal fields are always "external" and they are not even a subject for quantization. 
Therefore, we have to admit that the equations of the light front QCD evolution cannot have the first integral of 
the light-cone momentum. However, momentum conservation plays an important role in the derivation of the AP 
evolution equations, being solely responsible for the so-called +- prescription which regulates the collinear singularities 
of the splitting kernels. With momentum conservation the parton model is well motivated. Otherwise, we face the 
problem of identifying the subject of the QCD evolution itself. 



C. Evolution of classical fields and the structure function Fl 



To leading order, the standard evolution equations completely disregard the dynamics of the longitudinal fields 
in the QCD evo lution, thus missing a physically important part of the dynamical process. As it follows from the 
expression (2.37) for the longitudinal structure f uncti on Fl, the function Q is solely responsible for the scaling 
violation of Fl when the quarks are massless. Eq. (3.16) indicates that Q depends on p t only via the running coupling 
9r(Pt) ~ l/i°5(pf /A 2 )- Therefore, it is easy to estimate the <5 2 -dependence of the non-singlet longitudinal structure 
function Fl for the case of massless quarks, 



Q 2 F L 



4Q 2 



dp 2 



ln(p 2 /A 2 ) 



A 2 H(f) 



-A 2 E X ( 



ln (^)) 



4Q 2 



ln(4Q 2 /A 2 ) 



Q 2 a s (Q 2 



(3.29) 



where li(x) is the integral logarithm. Therefore, Fl is not strongly suppressed (i.e., by the powers of A 2 /Q 2 ) with 
respect to F 2 . In the standard OPE approach, the non-singlet Fl appears only in the next-to-leading order. Indeed, 
if the longitudinal fields were eliminated from the initial data by putting the partons on mass shell in a perturbative 
vacuum, at least one extra emission is needed to regenerate the longitudinal component at the level of the coefficient 
functions. For this reason, the OPE-based result contains additional <5 2 -dependence due to the scaling violation of 
the F2(x,Q 2 ) and G(x,Q 2 ). The OPE-based predictions for large x, when one might ex pect th e stro ngest effect from 
the longitudinal fields, systematically lie below the data. In our case, according to Eqs. (3.16)— (3.18), the geometry of 
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the static fields depends on the overall radiation above the currently probed value of the x f \ the entire range of Q 2 is 
integrated. An active involvement of the longitudinal (classical) fields in the QCD evolution exp lains the enhancement 
of the longitudinal structure function Fl due to the quark masses, the first two terms in Eq. ( 2.37 ); the heavier the 
quark, the more static is its field. 



IV. MASTER EQUATIONS OF QCD EVOLUTION 

Equations of motion of the quark and gluon fields do not allow one to eliminate longitudinal fields from the evolution 
equations by fiat. At most, one can use a renormalization procedure which absorbs their effects into the definition of 
some composite objects. The simplest objects of this kind are known as the particles in asymptotic states. Once this 
normalization is done, the only consistent way to proceed is to treat the entire problem as a problem of scattering. 
The initial state should be introduced as an asymptotic one. Thus, we are compelled to introduce the parton picture. 
The proton has to be presented as a system of free (valence?) quarks long before its interaction with the electron. 
Even the concept of "before" , which is so fruitful in the inclusive approach, becomes ambiguous. Instead of an initial 
picture of the inclusive e-p interaction, we must use the picture of the electron-parton scattering. In fact, we restrict 
ourselves to the study of the interaction between the electron and a single parton which is defined as an elementary 
excitation above the perturbative vacuum. The initial state of this parton has to be prepared. This is done implicitly 
by introducing a factorization scale. In fact, the factorization scale /ip undermines the basic idea of the inclusive 
measurement in e-p scattering. To comply with quantum mechanics, this parameter has to be explicitly measured. 
Moreover, this measured value has to be used in theoretical calculations! The recently discovered HERA events with 
rapidity gaps are good candidates to allow for this kind of measurements. 

This qualitative analysis promp ts fu rt her a pproximations. In order to obtain the master equations, 02 & n d W2 must 
be eliminated from the equations (3.9)— ( 3.13 ) without any discussion of the accuracy. After that, one can express <7q, 



a m , and W\ via the observables of DIS, and integrate the equations over p . We end up with two integral equations 
for quarks, 



( dq f (x,pj) 
\ dp\ 



9 2 (Pt) f'dy 
(2tt) 3 L y 

(k 2 + mj)(l + z 2 ) + 2m 2 f (l -z)\ d q f (y, k 2 t 



LL 
2 1 



[H(p,k lZ ) + {l- z)m 2 } 2 
{l-zf + z 2 f 1 

H(p, fc, z) + m 2 f 



dk 2 



d 2 k{ C F 



1 z 2 + l 

H(p, k, z) + (1 - z)m 2 1-2 
dM f (y,k 2 ) 



2m 2 z(l — z 



[H(p,k,: 



[H(p,k,z) + (1- z)m 2 ] 2 



z(l - z) k 2 



m 2 Z (1 — z) 



[H(p, k, z) 



,212 



dk 2 

dG(y,k 2 ) 



dk 2 



(4.1) 



dM f (x,p 2 ) L ... , r , 

LL 



dp 2 



9 2 (p 2 t) f'dy 

v 



(2nf 



and one equation for gluons 
'dG{x,p 2 )\ g 2 {p\ 



dp 



LL 



(2tt) 3 



d 2 k 



C F 



H(p, k, z) + (1 — z)m 2 



(1-z) 



dqf(y,k 2 
dk 2 



+ z 



dM f (y,k 2 y 



dk 2 



dG(y,k 2 ) 



dk 2 



H (p, k, z) + zm 2 [H(p, k, z) + zm 2 



212 



+N C 



2z(l — z)m 2 f 
[H(p, k, z) + zm 2 ] 2 
1-z 



dk 2 



dk 2 



1 - z 



+ z(l-z) 



1 



[H(p, k,z) + zm 2 ] 2 

(1 - z)k 2 1 dG(y, k 2 ) 



H(p,k,z) H 2 (p,k,z) 



dk 2 



2(H(p,k,z) + m 2 ) 
z(l - z)(k 2 + m 2 ) \ (l-z) 2 + l 



dq f (y,k 2 ) 2z 2 {\-z)m 2 dM f (y,k 2 ) 



(4.2) 



(4.3) 



where H (j>, k, z) = z(p— k) 2 + (1 — z)p\ . The subscript LL stands for the Leading Logarithms, the destination point 
of our study. It indicates that the rate of the evolution is considered with the switched off static fields. The further 
simplification of these equations will result in the DGLAP equations. One can easily see that Eq. (4.2) describes 
evolution of the dynamical quark mass (in fact, M. is the imaginary part of the "pole mass" which is defined by 
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the retarded self-energy of the quark). Negative signs on the right side of Eq. (4.2) indicate that evolution of the 
effective mass leads to a faster decrease, the higher the Q 2 which is probed. This trend clearly supports our intuitive 
unde rsta nding of the interplay between the quark mass and the transverse momentum transfer. By examination, 
Eq. (4.2) is of the same type as the other two evolution equations and the solution for M.f must be the standard 
logarithmic exponent. There is no reason to neglect A4 f(x,p 2 ), unless m 2 ? <C p 2 . 

Equations (|4.l| )-([4~g|) are still more complicated than the DGLAP equations. Unlike the DGLAP equations, they 
are the integral equations for the rates of evolution and exhibit non-locality in transverse directions. The kernels of 
these equations are singular; except for the familiar spurious poles, 2 = 1, of the splitting functions, the nonlocal part 
of the kernels contains overlapping infrared singularities at z = 1 and k — p. As discussed above, this behavior is quite 
consistent with the crude picture of the infinite momentum frame. Its smoothing requires a physical cut-off, which 
can be natural only in a more gentle picture. The cut-off brings into play new large logarithms, like log(s/A 2 ), which 
are alien to the dynamics of the infinite momentum frame and cannot really be estimated by its internal means. Even 
xp, the main variable of the parton model, has to be sacrificed in order to incorporate the energy of the collision as 
a physical parameter of the theory. In fact, this parameter is vitally needed in order to bring into the theory the 
Lorentz contraction of the colliding hadrons as the measure of the initial resolution in the longitudinal direction. 

To eliminate the need for the infrared cut-off, once again, by ordinance, one should follow the strategy of hunting 
the leading logarithms, either log(Q 2 / A 2 ), or log(l/x). As it was shown in the previous section, it is inconsistent to 
address the \ow-x behavior without explicit account of the longitudinal fields (however, they are already neglected in 
Eqs. (4.1 )— (4.3) ). In order to keep track of the leading logarithms log(Q 2 / A 2 ), one should find a solution to the 
scattering problem in the Born's approximation and exponentiate it. This is known also as ordering by angles or 
ordering of the emission by the transverse momenta. To obtain the same result from Eqs. ( f4.1| )-(f4~3|), it is necessary 
to eliminate non-local effects in their kernels (by neglecting k t with respect to p t ) and to introduce the upper limit p\ 
for the integration dk 2 . The integration itself becomes straightforward and yields, 
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(4.6) 



Thus, the angular ordering results in drastic simplification of the evolution equations. Only the splitting kernels 
P qq and P gg are singular now. However, we cannot regularize these singularities by means of the (+)-prescription as is 
done in AP equations. Indeed, in the AP approach the (+)-prescription maintains two conservation laws, conservation 
of the light cone charge, j + , and conservation of the light-cone component of the momentum, p + , in the process 
of sequential splitting of the parton. The parton is supposed to be free and no interaction terms are included into 
the operators j + and p + . Though the dynamical mass M. / of the fermion does not enter these operators, they are 
independent partners of the parton densities q and G in the evolution equations. Unless m/ = 0, the conservation laws 
do not follow from the (+)-prescription. This fact has a very clear physical interpretation. Though the wave equation 
for the massive quark is of hyperbolic type and allows for the propagation of signals along the light-cone characteristics, 
these signals cannot correspond to the on-mass-shell states of heavy quarks. Only dynamical polarization waves which 
include massive fermion fields can propagate at speed of light. Even though we can include the effect of quark masses 
into the gluon propagator and, eventually, into the running coupling, this will lead to an excess of accuracy if the 
dynamical mass term Mf is not included into the evolution equations. 

Our result for the running coupling |20] coincide with the one obtained by Dokshitser and Shirkov flit] , 



g 2 ^ 2 ) g 2 (M 2 ) 16tt 2 



ft 



In- 



M 2 



In- 



M 2 



/(4: 

m 2 



M 2 
/(— 
/ 



(4.7) 
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where the function / is obtained from the one-loop gluon self-energy with the subtraction at an arbitrary space-like 
momentum p 2 < 0, 




and has the following limiting behavior in different domains of the transverse momentum, 

i4i» i ; /(4)^ in y4' i4i « i- ( 4 - 9 ) 

Ttl p z Tfl Tfl 2m 1 Tfl 

These formulae smoothly interpolate running coupling between the domains with different transverse momenta and 
provide different numbers of "active flavors" . However, it is important to emphasize that with the quark masses 
retained, the dynamical quark mass should be kept in the evolution equations and we still face problem of the infrared 
regularization of the evolution equations. 



V. SUMMARY 



The method of QFK allows one to extend the standard definition of the evolution equations beyond the scope 
of the parton model. The extended equations include new elements, the static (longitudinal) gluon field and the 
static (non-propagating) quark fields. The quark masses are accounted for also and the dynamical mass term in the 
evolution equations is sensitive to the change of scale in the same way as the quark and gluon structure functions. 
Though the role of the new terms is perfectly understood physically, the extended set of equations suffers from the 
infrared-singular terms which cannot be regulated in the the standard scheme where the proton is analyzed in the 
infinite momentum frame. 
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